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Abstract. We report the discovery of the new dipping X-ray source, XMMU J004308.6+411247, in M31, during a systematic 
search for periodicities in XMM-A'ewton archival observations. During the 2002 January 6 observation, the dips recur with 
a 107 min period and the source count rate is consistent with zero at the dip minimum. Dips with the same modulation and 
period are also observed during the XMM-Newton observations carried out on 2000 June 25, 2001 June 29 and the Chandra 
observation of 2001 October 5. The dips of XMMU J004308. 6+41 1247 show no evidence of energy dependence. The average 
X-ray flux of XMMU J004308. 6+41 1247 is nearly constant across dilferent observations (~ 10^' erg s"' for an assumed M31 
distance of 780 kpc in the 0.3-10 keV band); the spectrum is well fit by an absorbed power law with a photon index ~ 0.8 or 
an absorbed Comptonization model. The photo-electric absorption is consistent with the Galactic value in the source direction. 
If XMMU J004308. 6+41 1247 is located in M31 its properties are consistent with those of dipping low mass X-ray binaries in 
the Galaxy. Present observations do not allow to distinguish between dips and eclipses. 

The possibility that XMMU J004308. 6+41 1247 is a foreground X-ray source cannot be ruled out at present; in this case 
the source might be a magnetic cataclysmic variable. 
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1. Introduction 
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Q ' Being the closest large spiral (Sb) galaxy to our own, the Andromeda Galaxy (M31) provides us with a prime opportunity to 
study X-ray emission in a Galaxy similar to the Milky Way both in morphology and size. The sources in M31 are observed 
^ at the nearly uniform and well known distance of ~780kpc (Stanek & Garnavich 1998, Macri et al. I2001t . Moreover, owing 
to its inclination angle of ~ 77° and high galactic latitude (ZII= 121.174312, WI= -21.573022) sources are viewed through a 
K*" substantially lower absorption column (A^h ~ 7 x 10^" cm^. Dickey & Lockman 1990") than for sources in the Galaxy. Because 
of this moderate and fairly uniform extinction, M3 1 sources can be studied and compared over an extended low energy band. 
5J] ■ Moreover, different stellar populations in the bulge, spiral arms and halo can be easily distinguished based on their location. 

The XMM-Newton mission ( Jansen et al. 1200 It yields an unprecedented opportunity to study the time variability and spectral 
properties of individual X-ray sources in nearby galaxies with increased throughput. Previous observations of M3 1 performed 
with Einstein (Trinchieri & Fabbiano ll991l van Speybroeck et al. \\919\ and ROSAT {Vhmm\, Forman & Jones fT993b revealed 
many point-like X-ray sources in its core (r < 5' ~ 1.1 kpc), but were too insensitive for their precise characterization. The two 
deep and extensive ROSAT-PSPC surveys of the entire M31 disk (Supper et al. 2001 ; Supper et al. , 199 7 \. though effective in 
mapping the Globular Cluster (GC) and Supernova Remnant (SNR) populations, proved inadequate for a detailed study of point 
sources in the core of the galaxy because of the limited spatial resolution. Spatial resolution was the main limitation also in the 
study of the bulge of M31 caiTied out with Beppo-SAX, the first to image M31 up to energies of ~10keV (Trinchieri et al. \199^ . 



The observations of the central part of M3 1 with the Chandra X-Ray Observatory (Wei sskopf ll988l l resolved the nuclear 
source seen with the Einstein and ROSAT-HRI into five point sources (Garcia et al. '2000) and detected and localized bulge point 
sources down to luminosities of ~ 2 x 10^^ erg s"' in the 0.3-7 keV band, with positional eiTors of less than 1" (Kong et a/. 120021 



Di Stefano et al. 120021 Kaaret E002l l. None of these observations, however, was long enough for a detailed timing analysis. The 
large collecting area and bandpass of XMM-A^ewfon afforded a more accurate characterization of both the global properties of 
the X-ray emission from M31 (Shirey et al. 2001 ) and individual X-ray sources in its bulge (Osborne et al. 2001 Barnard, Kolb 
& Osborne [20021 Trudolyubov et al. [2002 ). M3 1 was selected as an XMM-A^ewfon Performance Verification target (Jansen et al. 
1200 H and subsequently observed within the Guaranteed Time Program. The XMM-Newton observations of M3 1 (see section|5J 
are among the M3 1 X-ray observations with the longest total exposure; they are eminently suited for timing and spectral analysis. 
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Table 1. XMM-Newton observations of M31 used in this analysis. 



Observation 


Date 


^start 


^stop 


Exposure time MOS" 


Exposure time PN" 




(UT) 


(UT) 


(UT) 


(ks) 


(ks) 


1 


25/06/2000 


08:12:41 


20:59:33 


34.8 (30.1) 


31.0 (26.5) 


2 


28/12/2000 


00:01:33 


03:43:23 


12.2(12.2) 


9.9 (9.9) 


3 


29/06/2001 


06:15:17 


22:20:20 


54.9 (31.0) 


52.4 (31.0) 


4 


06/01/2002 


18:00:56 


11:52:53* 


63.0 (63.0) 


60.1 (60.1) 



" Effective exposures used in our analysis. The number in parentheses is the exposure after exclusion of time intervals with high background 

levels. 
* The following day. 



We started a program of systematic search for periodicities with Fourier techniques among all bright ( 5 x 10^^ erg s ') point- 
like X-ray sources detected in the XMM-Newton archival observations of M3 1 . Our search strategy is aimed at detecting periodic 
signals over the widest possible period range (see sectionl^. 

In this paper we report the discovery of a 107 min dip-like modulation (nearly 100% amplitude) in the 2002 January 6 light 
curve of the source XMMU J004308. 6-1-41 1247. The same modulation is clearly seen also in other XMM-A^evvfon observations 
of M31 and in the ACIS-S Chandra observation of 2001 October 5 (see section BTTl . 

In the following sections we describe the set of data we analysed (section |3, summarise the periodicity search technique 
(section|3j and illustrate the properties of the dip we discovered in the light curve of XMMU J004308. 6-1-41 1247 (section|3. Our 
results are discussed in section|5l 

Throughout the paper we adopt a distance to M31 of 780 kpc (Stanek & Garnavich ll998l Macri et al. 1200 1> and a Galactic 
column density in the direction of M3 1 of 7 x 10^** cm"^ (Dickey & Lockman ll990> . 

2. Observations and data analysis 

We used the four archival XMM-Newton observations of the core of M3 1 listed in Table The 2000 June 25 observation 
(observation 1 in Tabled was part of the Performance Verification Program (Jansen et al. 2001 Shirey et al. 120011 Osborne 
et al. I2001> . The other three observations (observations 2,3,4) were carried out on 2000 December 28, 2001 June 29 and 2002 
January 6 respectively as a part of Guaranteed Time Program. We used data from the three European Photon Imaging Camera 
(EPIC) instruments at the focus of the three co-aligned 1500 cm^ X-ray telescopes on board XMM-A^ewfo«: two EPIC MOS 
detectors (Turner et al. 2001 1 and one EPIC PN detector (Struder et al. 2001 1. Each MOS detector is sensitive in the 0.2-lOkeV 
band, while the PN detector is sensitive in the 0.1-15keV band. In all observations, the three quoted detectors were operated 
in the full window mode (~ 30' diameter field of view) with medium (observations 1, 2 and 3) and thin (observation 4) optical 
blocking filters. We reduced EPIC data with the XMM-A^ewton Science Analysis System (SAS) version 5.4.1, by performing 
standard screening of the EPIC data in order to exclude time intervals with high background levels and/or solar flares. This 
resulted in shortened exposures of observations 1 and 3 as indicated in Table We applied the automatic source detection and 
analysis pipeline summarised in the following section to all the four observations. Observations 1, 3 and 4 have high enough 
statistics to carry out a detailed timing analysis of more than half of the individual point sources in them. Observation 2 has a 
somewhat lower exposure. Thus, section0]present only the results we obtained from observations 1, 3 and 4. In Fig.^the central 
part of the MOS 1 image of M3 1 extracted from the 2002 January 6 observation is shown. 

3. Searcli strategy 

For each observation listed in table[l]we filtered PN, MOSl and MOS2 events in the energy band 0.2-15 keV, and merged the 
events of the two MOS detectors into a single MOS1+MOS2 event list (task merge). We then produced images of the PN and 
MOS1H-MOS2 fields binned in square pixels of 2" x 2". Note that this is smaller than the FWHM of the point spread function in 
the center of the field for both PN and MOS, which is ~ 6" for the PN and ~ 4" for the MOS (XMM-Newton Users' Handbook 
V2. 1 1200 ^. We detected all sources above the 5cr threshold in the two images (task ewavelet). For a description of the wavelet 
based source detection algorithm used by the task ewavelet see Damiani et al. 199? For every detected source this task yields 
the source position and extension, corresponding to the standard deviation of an equivalent Gaussian luminosity profile, i.e. the 
radius containing ~66% of the source counts. In the longest exposure observation (observation 4) we found ~ 350 sources in the 
PN image and ~ 180 in the MOS image. We removed from the source lists some spurious sources close to the image edge, CCD 
borders and bad pixel lines (this was important especially in PN images). 

After that we extracted light curves of all detected point sources and applied the barycentric correction to the photon arrival 
times. We used an extraction radius of two times the source extension calculated by ewavelet (i.e. 2cr of the equivalent Gaussian 
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Fig. 1. Image of M31 from the XMM-A^ewfon observation carried out on 2002 January 6 (MOSl data). The source 
XMMU J004308.6+41 1247 is marked by the extraction circle that we used, containing 98% of the source counts. 



luminosity profile; this contains about 98% of the source counts) whenever there was no superposition with nearby sources, and 
an extraction radius of one time the source extension in all other cases. 

Finally we considered all extracted light curves containing more than 200 photons (including background) in order to have 
enough statistics for timing analysis. In this way, for instance, we reduced the total number of analysed sources to ~140 (PN) 
and ~90 (MOS) in observation 4. Detailed discussion on the complete set of sources that we analysed will be presented in a 
forthcoming paper. 

The light curves obtained were searched for pulsations by using the power spectrum technique developed by Israel & Stella 
J1996> . In order to maximize signal detection sensitivity we calculated a single power spectrum with the highest Fourier resolution 
Av = 1/Ar (with Ar the time span covered by the observation). In each power spectrum we searched for peaks above the 
detection threshold up to the Nyquist frequency VNyq = 1/2 Af, with the binning time Af equal to the intrinsic time resolution of 
the EPIC cameras (i.e. 73.4 ms for the PN and 2.6 s for the MOS). The detection threshold for peaks in a single power spectrum 
containing Ntriai - ISJ jlAt frequencies was set such that the probability of exceeding it by chance in a total number of trials 
equal to Nmai times the number n of light curves analysed in the image was e = « * Nmai * esingie trial, with e - 2.326 x 10"'* 
(corresponding to 3.5cr in a Gaussian approximation). Thus, the detection threshold for the analysis of a single power spectrum 
was given by e/n. In this way, the minimum detectable amplitude (or pulsed fraction) for a sinusoidal signal (i.e. the sensitivity 
of the search to coherent periodic signals) for each source can be derived according to the prescription given by Vaughan et al. 



4. Results 

4.1. Timing analysis 

The systematic search described in section |3] revealed both the 865 s pulsation of the transient super-soft source 
XMMU J004319.4+41 1759 (Osborne et al. 1200 U visible in the 2000 June 25 observation (observation 1) and quiescent in 
all other archival XMM-A^evv'fo« observations, and the 2.78 hr dips in the 2002 January 6 observation (observation 4) of 
XMMU J004314.1+410724in the globular cluster Bo 158 (Trudolyubov et a/. .2002j . 
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Fig. 2. Power spectrum of the combined PN and MOS light curves of the source XMMU J004308. 6+41 1247 corresponding to 
the 2002 January 6 observation. The power spectrum is normalized according to Leahy et al. il983i 



We also found a significant peak at 1.525879 x 10"^ Hz in the power spectrum of XMMU J004308. 6+41 1247 extracted 
from observation 4 MOS data. The source position we obtained, RA = 00H3'"08';61 1 DEC = 4ri2'47'.'58 (equinox 2000), is 
to be considered accurate to better than 4". The power spectrum peak has a significance of ~ 4.3cr after taking into account the 
number of trial frequencies for this source, as well as for the other sources in the MOS image of observation 4 for which the 
periodicity search could be carried out. Visual inspection of the power spectra of the 2002 January 6 MOS and PN light curves 
of XMMU J004308. 6+41 1247 clearly reveals the presence of the second and the third harmonics of the signal. As a further 
confirmation, we calculated the power spectrum of the combined PN and MOS light curves (using the XRONOS package, 
version 5.19) and obtained a significant increase in the power of the signal and its harmonics. In Fig.|2the power spectrum of the 
PN+MOS light curve of XMMU J004308.6+41 1247 is shown. The signal fundamental plus two harmonics can be easily seen. 
An epoch folding search was carried out on the PN, MOS and PN+MOS light curves of XMMU J004308.6+41 1247 from 2002 
January 6 observation aiming at a more accurate determination of the period value. Following the prescription of Leahy ( 1987 ) we 
obtained an average best period of 6420 ± 80 s (error estimate at 67% confidence level). For the background count rate estimation 
we used an annular region around the extraction circle of XMMU J004308. 6+41 1247 (see also section The background 
subtracted PN+MOS light curve folded at a period of 6420 s is shown in Fig.|3(panel a). The reference epoch Tq corresponding 
to the minimum of the modulation is given in Table|2] together with the source count rate. The minima are consistent with a zero 
count rate. The dips last ~25% of each cycle and their fall and rise are fairly continuous. A fit of the folded light curve in Fig|3l 
with a simple model consisting of a constant and a Gaussian with negative normalization centered at the minimum gives a nearly 
100% amplitude of the dip with respect to the out-of-dip intensity and a cr = (8.3 ± 1.2) x 10"^ (phase units). See Table|2 The 
same results are obtained from the combined PN+MOS light curve restricted to the 0.2-1 keV and 1-15 keV energy bands and 
folded at the 6420 s period (see Fig.|3 panels b and c). The fit to the folded light curve gives a dip amplitude of 104%±27% 
and 95%+ 12% in the 0.2-1 keV and in the 1-15 keV energy band respectively, fully compatible with the 98%±1 1% amplitude 
obtained in the 0.2-15 keV energy band. We thus conclude that the data provide no evidence for an energy dependence of the 
modulation. 

We searched for a modulation at a period around 6420 s also in all the other HMM-Newton observations of 
XMMU J004308. 6+41 1247, using combined PN+MOS light curves of the source in all cases in order to improve statistics. 
An epoch folding search yields a best period of 6400+440 s in observation 1 and 6425+200 s in observation 3. The errors are a 
factor 3-5 larger than in observation 4 and the best period values are consistent with one another. Therefore we are justified in 
adopting the same period of 6420 s to fold all the light curves. We folded the XMMU J004308.6+41 1247 light curves exti-acted 
from XMM-Newton observations 1 and 3 at the 6420 s period and fitted the folded light curve with the same model used for 
observation 4 (a constant and a Gaussian with negative normalization). We confirmed the presence of the signal at the same 
period and with a modulation amplitude compatible with that previously found. Results are shown in Fig.|3(panels a and b) and 
summarised in Tabled 
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Fig. 3. Panel (a) shows the combined PN+MOS hght curve of XMMU J004308.6+41 1247 from the XMM-Newton observation 
of 2002 January 6 (observation 4), in the energy band 0.2-15 keV, and folded at the best period of 6420 s. The solid line represents 
the fit of the data with the model described in the text (a constant minus a Gaussian). Panels (b) and (c) show the results obtained 
by folding at 6420 s the combined PN+MOS light curve in the 0.2-1 keV and in the 1-15 keV energy band, respectively. Each 
light curve is background subtracted and normalized to the average source count rate value given in the top right corner of the 
panel. 



In order to further investigate XMMU J004308.6+41 1247 we used archival Chandra and TJO^Ar observations of M31. We 
selected the archival Chandra observations with the longest exposure, i.e. that carried out on 2001 October 5 with the Imaging 
Spectrometer ACIS-S (Garmire et al. L992, Bautz 1998 and references therein) for an eflFective exposure of ~37.7 ks and the 
2001 November 1 observation with the HRC-I in imaging mode for an effective exposure of ~46.8 ks (Murray et al. I1997t . 
We used standard event processing and filtering procedures in the CIAO package. Version 2.2.1, and extracted light curves of 
XMMU J004308.6+41 1247 from circular regions of diameter about twice the FWHM of the local Point Spread Function (PSF), 
in the 0.3-lOkeV band. We folded the XMMU J004308.6+41 1247 hght curves obtained from Chandra data at the 6420 s period. 
The dips are clearly detected in the ACIS-S Chandra observation carried out on 2001 October 5 (see panel c in Fig.|4j and 
their properties are consistent with those of the XMM-A^evvfon observations (see Table |2j. The best period obtained with an 
epoch folding search on the ACIS-S Chandra observation of 2001 October 5 is 6437±200s, consistent with the period during 
observation 4. There is marginal evidence for the 107min signal also in HRC-I Chandra observation of 2001 November 1 (see 
Fig-E} though the epoch folding search do not show any significant peak. 

We then considered the longest exposure 7?(9&4r archival observations of M3 1 containing XMMU J004308. 6+41 1247 within 
ten arc-minutes from the center of the field of view. These are the 1996 January 7 HRI observation (85.5 ks exposure) and the 
1991 July 27 PSPC observation (30 ks exposure). The ROSATlight curves of XMMU J004308.6+41 1247 obtained contain too 
few photons (~ 60 and ~ 80 net source counts respectively) to confirm the presence of the dips. 

Note that the four folded light curves of XMMU J004308. 6+41 1247 obtained from the XMM-Newton and Chandra 
data can be reasonably well fitted also by a square-wave. For this we adopted a model of the form Y{F) - 
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Fig. 4. We show the combined PN+MOS light curves of XMMU J004308.6+41 1247 in the band 0.2-15 keV folded at 107 min 
period extracted from XMM-Newton observations of 2000 June 25 (a) and 2001 June 29 (b). The fit with a simple model 
consisting of a constant and a Gaussian with negative normaUzation is also shown (solid hne). In panel (c) we show the folded 
light curve of an ACIS-S Chandra observation of the source of 2001 October 5. Each light curve is background subtracted and 
normalized to the source average count rate. This is given in the top right comer of each panel. 



Table 2. Parameters of the observed modulation 





Mission 


Date (UT) 


To" (MJD) 


SOURCE* (counts/s) 
PN MOS 




A" 


1 

3 
4 


XMM 
XMM 
XMM 


25/06/2000 
29/06/2001 
06/01/2002 


51720.5201+1.5 x lO'^ 
52089.4458+7 x 10"" 
52280.7957±9 x 10"" 


1.7+0.1 0.98+0.08 
1.4+0.2 1.0+0.1 
1.32+0.09 1.08+0.06 


0.07+0.015 
0.05+0.015 
0.08+0.01 


88%+20% 
115%+18% 
98%+ll% 




Mission 


Date (UT) 


To "(MJD) 


SOURCE* (counts/s) 
ACIS 


a' 


A" 




Chandra 


05/10/2001 


52187.057+2 x 10"' 


0.99+0.05 


0.07+0.02 


92%+ 17% 



" Reference epoch corresponding to the minimum of the modulation 

* Source count rate in units of 10"^ counts/s in the 0.2-15 keV band for the XMM-Newton observations and in the 0.3-10 keV band for the 
Chandra observation. The MOS colvuim refers to the MOS1+MOS2 data. 
Standard deviation of the Gaussian dip in phase units. Errors are at the 90% confidence level. 
Amplitude of the modulation with respect to the out-of-dip intensity. Errors are at the 90% confidence level. 



D [arctan(A(Fi - F))/7: + arctan(A(F - F2))/7t] + C, where A, Fi, F2, D and C are the parameters of the model and F repre- 
sent the phase. For very large values of A(F2 - Fi) (> 1000) this expression produces an echpse-type behaviour with sharp 
ingress and egress at phases F\ and F2, respectively. When A(F2 - Fi) is smaller the ingress and egress become progressively 
smoother, Fi and F2 being the half depth phases. When A(F2 - Fi) i 100 the profile turns into a dip-like one with no flat 
minimum phase interval. 
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Fig. 5. Background subtracted and folded light curve of XMMU J004308. 6+411247 from the Chandra HRC-I observation of 
2001 November 1 . The average source count rate, is given in top right corner of the panel. 



We fit all our folded light curves to the square-wave model with A frozen to 10'*, the being insensitive to changes in A 
within one order of magnitude from this value. The results are shown in Fig. |6l The square-wave model fit provides a reduced 

slightly closer to 1 than the Gaussian fit, but uncertainties in the parameters D and C (needed for the eclipse depth estimation) 
are large. According the square-wave model fit the XMM-Newton 2002 January 6 light curve profile shows an eclipse lasting 
0.19+0.03 phase units with a 92%±38% decrease of the source flux during the eclipse (90% confidence level errors). In the 
XMM-Newton 2000 June 25 light curve the eclipse lasts 0.21+0.06 phase units and produces a 84%+45% flux decrease; in the 
XMM-Newton 2001 June 29 it lasts 0.125+0.005 and gives a 104%+60% flux decrease and finally, in the Chandra 2001 October 
5 observation the eclipse lasts 0.21+0.03 and gives an 89%±49% flux decrease. 

Because of these results we conclude that poor statistics does not allow us to distinguish between a smooth profile of the 
folded light curve profile (as the one represented by the Gaussian model) or a sharp eclipse-like light curve profile. 



4.2. Spectral analysis 

We extracted phase-averaged spectra of XMMU J004308.6+41 1247 in the 0.2-15 keV energy band from the PN and MOS data 
of all the XMM-Newton observations in which the 107 min modulation was detected (observations 1, 3 and 4). We generated 
the relevant spectral response matrices (SAS tasks rmfgen and arfgen). The background spectrum was extracted from the same 
annular region around the source extraction circle used for the background count rate estimation. These spectra were analysed 
simultaneously by fitting the following spectral models: absorbed power law, absorbed power law and blackbody, absorbed power 
law with exponential cut off, absorbed thermal bremsstrahlung, absorbed Comptonization model (Titarchuk 1994 Titarchuk & 
Lvubarskii ll995> with spherical or disk geometry. Given the relatively poor statistics of our spectra, it is perhaps not surprising 
that a simple absorbed power law provides a statistically acceptable fit (see below). Nonetheless we fit also a larger sample of 
two components spectral models that are often used to describe the spectra of various classes of Galactic pulsating X-ray sources 
in order to try to draw some closer analogies. The spectral fitting was performed using the same model parameters for the PN, 
MOS 1 and M0S2 phase-averaged spectra. We allowed for independent normalizations of the three spectra in order to account 
for calibration differences between the EPIC cameras. Results from the spectral fitting are given in Table|3l Only models yielding 
a reduced less than 2 are reported. The model consisting of an absorbed power law with an exponential cut off was rejected 
because it yields a cutoff energy greater than 100 keV, i.e. a spectrum virtually indistinguishable from a power law in the XMM- 
Newton energy band. Note that the absorbed Comptonization model (Titarchuk . 1994) gives an acceptable fit only by assuming a 
disk geometry for the Comptonising region. 

The model consisting of a simple absorbed power law, gives the best reduced in all observations, but provides an upper 
limit to the absorption column at 95% confidence level that is lower than the Galactic value of 7 x 10^" cm"^ derived from the 
HI distribution map by Dickey & Lockman I 19901 The latter has a very limited angular resolution (1° x 1°). It was estimated that 
small scale structure in the HI distribution in the sky introduces errors in the column density estimates of ±1 x 10^" cm"^ (90% 
confidence level, see Elvis et al. 1986 Appendix B). The upper limits to A^h derived from the absorbed power-law fit are below 
the 90% confidence interval for the Galactic value. The fits with the Comptonization model has the same problem in two out 
of three observations. The inclusion of a black body component in the absorbed power law model (this is done in analogy with 
the X-ray spectra modeling used for several Galactic dipping sources, see e.g. XB 1916-053, Church et a/. 1 19971 XB 1254-690, 
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Fig. 6. Square-wave fit of the background subtracted and folded light curves of XMMU J004308.6+41 1247 from the XMM- 
Newton and Chandra observations showing the 107 min pulsation. The satellite, observation date and average source count rate 
are given in each panel. The data are the same as in Fig. Impanel (a) and in Fig. Impanels (a), (b) and (c). 



Smale, Church & Baluciriska-Church'2002' X 1624^9, Church & Baluciriska-Church'1995^ gives plausible upper limits to the 
absorption column, but is not required statistically since the power-law fit is akeady acceptable. The phase-averaged spectra of 
XMMU J004308.6-H41 1247 obtained from the 2002 January 6 observation are shown in Fig.|7l 

In order to investigate further the low-level absorption indicated by the above fits, we also fit groups of PN, MOS 1 and MOS2 
phase-averaged spectra with the three above-discussed models, while freezing A^h at the Galactic value (which is to be considered 
a lower limit to the absorption for a source located in M31). We obtain reduced Jt'^ values of 1.2 (power law), 1.03 (power law 
plus blackbody) and 1.1 (comptonization) for observation 1, of 1.4, 0.96 and 1.2 for observation 3 and of 1.3, 0.76 and 1.06 for 
observation 4 respectively and parameter values consistent with the corresponding ones listed in Table |3 to within the errors. 
Note also that the upper Umits to A^h increases considerably (up to values consistent with the Galactic value in some cases) 
when allowing for independent absorption columns in the PN, MOS 1 and MOS2 spectra. The 2002 January 6 PN, MOS 1 and 
MOS2 phase-averaged spectra fit in this way with the absorbed power-law model gives a^^ of 48.1 (47 d.o.f.), a photon index 
of 0.79+°-^5 and A^h < 4 x 10^" cm^^ for the PN spectrum, A^h < 6.8 x 10^" cm^^ for the MOSl spectrum, A^h < 14 x 10^" cm^^ 
for the MOS2 spectrum. Fitting again the three spectra with the absorbed Comptonization model (disk geometry), we obtain a 
of 37.20 (45 d.o.f.), an input photon temperature of kT^ < 0.12keV, Comptonising electron temperature of kT^ - 3.3^^ ^ keV, 
and optical depth t - 12.5;^] and A^h < 6.8 x 10^" cm"^ for the PN spectrum, A^h < 1 1 -5 x 10^" cm"^ for the MOS 1 spectrum, 
A^H < 1.7 X 10^*' cm"^ for the MOS2 spectrum. Based on this analysis we conclude that the low 95% upper limits to A^h obtained 
may be related to calibration uncertainties of the EPIC cameras. It is worth noting that XMMU J004308. 6-1-41 1247 is not the only 
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Parameter 


Observation Date 






2000 June 25 2001 June 29 


2002 January 6 




(observation 1) (observation 3) 


(observation 4) 


Absorbed Power Law [wabs*powerlaw] 



Photon Index 
Flux 
xHd.o.f) 



< 0.058 
0.84!»1^ 
1 71+"" 
21.53(20) 



(i-76!S:;^) 



< 0.023 

75+'' '3 

1 56+" " 
27.48(25) 



(1-56^;^) 



< 0.012 

O+0.07 
^-0.07 



l.( 



48.13(49) 



N, 



H 



Photon Index 

Flux ' 
Flux (black body) 
Flux (power law) ' 
xHd.o.f) 



(1.68--) 



Absorbed Power Law plus Black Body [wabs*(powerlaw+blackbody)] 



<0.32 




0.07 " 




< 0.098 




1.81+?^ 




2 14+" 




1 7Q+0.82 




2 3+^-'' 




9 07+0.96 
^■"'-04 




2 55+' "' 




1 69+°'*' 1 


f 1 Q 1 +0.901 


1 1 50+'''''' 1 


( 1 CQ+0.49"! 


1 1 7Q+0.37 1 

' -^-0.23 1 


[l.79: 


1 14+0-63 1 

^•^ -0.29 < 


1 1 1 c+0.64' 
1-^-0.29; 


1 1 20+""" 1 


' 1 20+""^ 1 
, i.Z,U„ 1,1 


1 1 21+°-'"' 1 


1.21: 


55+"-22 1 




1 30+" °* 1 


'0 36+" °^ 1 

(,"■-'"-0.06) 


1 52+''-'''' 1 

' "--^^-0.04 1 


[o.58: 


17.25(16) 




21.09(22) 




34.77(45) 





+0.37 \ 

-0.23 / 
+0.32 \ 
-0.19/ 
+0.05 \ 
-0.05 / 



Absorbed Comptonization Model [wabs*comptt] 



kT, f 

T = 

Flux = 
xHd.o.f) 



< 0.096 
<0.23 

T 9+49.2 

1 1 5+^-5 

1 7Q+0.72 
^•'-'-0.45 

20.20(18) 



< 0.044 

< 0.17 

2 84+2--^ 

^•°^-0.85 

14 1+'-' 
1 60+" " 
22.81(23) 



(l-60!-l) 



< 0.037 

< 0.15 

3 2+--2 

J3+2.0 

1 77+0.56 
'-0.23 

40.6(47) 



(1.78!°-) 



^ Results of simultaneous fit of EPIC-PN, MOSl and M0S2 phase-averaged spectra in the 0.2-15 keV energy band. Parameter errors 

correspond to 90% confidence level. 
^ Upper limit at 95% confidence level to the equivalent absorbing hydrogen column density in units of 10"^ cm"^ 

Observed flux in the 0.3-lOkeV energy range in units of lO"''' ergs"' cm"^. Average value from PN, MOSl and M0S2. The value of the 

unabsorbed flux is in parentheses. Errors are calculated from the uncertainty in the normalization at the 67% confidence level. 

Blackbody temperature in keV 

Upper limit to the temperature of soft photons in keV 
^ Electron temperature in keV 
^ Thomson optical depth for the disk geometry 
^ Parameter frozen to the Galactic value 



source we found showing this anomalous behaviour. Also the bright source XMMU J004232.0+41 1314show spectra with typical 
photon index ~ 1 and A^h < 0.6 x 10^" cm^ at 95% confidence level in the 2002 January 6 observation. In this case the upper 
limit to A^H grows to 0.9 x 10^" cm^ and to 1.2 x 10^° in MOSl and MOS2 spectra respectively when allowing for independent 
A^H parameters in the fit (we did not extract PN spectra because in the PN image the source is more than 50% covered by CCD 
edges and hot pixel lines). XMMU J004232. 0+41 1314 corresponds to source r2-32 in the ACIS-I Chandra catalogue by Kong 
et al. (ED02i. In the quoted paper (Table 8) the source was listed among the 20 brightest sources in M31 and its spectrum was 
well represented by a power-law spectral fit in the 0.5-lOkeV band with photon index 1 .3+qj and A^h < 5 x 10^** cm^ at the 95% 
confidence level, i.e. an upper limit to the absorption column lower than the Galactic value. 

Based on the absorbed flux values in the 0.3-lOkeV band given in Table |3l we can conclude that the luminosity of 
XMMU J004308. 6+41 1247 was roughly constant over the ~ 1.5 yr time span covered by XMM-Newton observations, at a 
level of ~ 1 X 10" (d/dMnf erg s"', where d is the actual distance of XMMU J004308.6+41 1247 fr om us and duii the M31 
distance. 

In order to investigate changes in the energy spectrum during the dips, we divided the PN, MOS 1 and MOS2 events into 
two phase intervals corresponding to the out-of-dip (from phase 0.1 to phase 0.85) and dip intensity (from phase 0.85 to 0.1). 
We extracted spectra from events in the two quoted phase intervals, thus obtaining out-of-dip and dip spectra. The MOS 1 and 
MOS2 dip spectra of observations 1 and 3 had less than 30 counts each (including background) and therefore were not used in 
any spectral fit. For each observation we fit the PN out-of-dip and dip spectra simultaneously with the phase-averaged spectrum, 
using two different models: an absorbed power law and an absorbed Comptonization model. Independent normalizations and 
absorption columns for the three spectra were used, while the other parameters were frozen to the best fit values in Table|3] To 
within the obvious statistical limitations this approach is justified by the fact that simultaneous fits of the PN, MOS 1 and MOS2 
out-of-dip spectra alone give results consistent with those obtained for the phase-averaged spectra listed in Table |3] A similar 
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Table 4. Best-fit model parameters of energy specti-a of XMMU J004308.6+41 1247 " 



Parameter 






Observation Date 








2000 June 25 


2001 June 29 


2002 January 6 




(observation 1) 


(observation 3) 


(observation 4) 




PN data 


PN data 


PN,M0S1 and M0S2 data 


Absorbed Power Law [wabs*powerlaw] 


Nh (out-of-dip) 


< 0.023 




< 0.037 




<0.01 




NH(dip) " 


< 9.7 




<0.93 




< 0.22 




Flux (out-of-dip) ° 






1 i.89t;^ 1 






(2.02!°:gn 


Flux ( dip) ' 


i.orjJJl 1 


(l.09!°]^j 






0.69!«;;5 1 


(0.69t«) 


^2(d.o.f) 


39.74(26) 




27.48(24) 




85.85(95) 




Absorbed Comptonization Model [wabs*comptt] 


Nh (out-of-dip) 


< 0.044 




< 0.048 




<0.03 




NH(dip) " 


< 32 




<0.90 




<0.19 




Flux (out-of-dip) 


2 14+"-'''* 1 


;2.l4!!!;t^) 


1 1.93!°;^^ 1 




2.11!«;]^ 1 


!2.13t?^) 


Flux ( dip) 


1 03+°'*'* 1 




1 l.ll!i;g 1 


^1.11 029 j 




(0.73!-0 


^'(d.o.fl 


33.12(25) 




22.81(23) 




70.79(94) 





" Results of simultaneous fit of the out-of-dip, dip and phase-averaged spectra of source XMMU J004308. 6-1-41 1247. Spectral parameters 
are frozen to the values indicated in Table|3| 

Upper limit at 95% confidence level to the equivalent absorbing hydrogen column density in units of 10~^ cm"^ 

Model observed flux in the 0.3-10 keV energy range in units of 10"'^^ ergs"' cm"-. The value of the unabsorbed flux is in parentheses. 
Errors are calculated from the uncertainty in the normalization at the 67% confidence level. 



check cannot be done with dip spectra alone because of poor statistics, and the usage of the same spectral model fitting out-of-dip 
and phase-averaged spectra represents the simplest assumption. Results are shown in Table |3 Note that for the 2002 January 6 
observation also the MOS 1 and MOS2 phase-averaged, out-of-dip and dip spectra were fit simultaneously to the PN data. In this 
case we used nine independent normalizations for the nine spectra and three independent A^h parameters, one for the three (PN, 
MOS 1 ans MOS2) phase-averaged spectra, one for the out-of-dip spectra and one for the dip spectra. The PN dip and out-of-dip 
spectra from the 2002 January 6 observation are shown in fig.|8ltogether with their best fit power law model. The resulting 95% 
confidence upper limit on the average increase in the (neutral) absorbing column density A^h during the dips was 9.7 x 10^^ cm~^ 
in observation 1, 9 x 10^' cm"^ in observation 3 and 2 x 10^' cm~^ in observation 4. (The tighter limit obtained for observation 
4 clearly reflects the longer exposure time in that observation). However the above limits should be taken as indicative, since in 
order to increase the S/N of the dip spectrum was integrated over the whole duration of the dips (as opposed to a shorter phase 
interval around the dip minimum). 



4.3. Other Observations 

An X-ray source with a position consistent with that of XMMU J004308. 6-1-41 1247 was detected with the ROSAT-HRI (source 
[PFJ93]74, Primini, Forman & Jones 1993 1, ROSAT-PSPC (sources 2RXP J004308.5+41 1251 and 2RXP J004308.3-I-41 1249 in 
the Second ROSAT PSPC Catalogue; source RX J0043.1+4112 in the M31 second ROSAT PSPC survey Supper et al. 2QW^, 
and Chandra (source CXOM31 J004308.4+41 1247 in Kaaret 2002 source CX0M31 J004308.6+41 1248 in Kong et al. .20021 
source CXOM31 J004308.6-i-411250in Wifliams et al. ^B^ . 

The observed flux of XMMU J004308. 6-1-41 1247 from the latter experiments is roughly constant over a 10 years time span. 
^From the published values of the source luminosity we derive values of the integrated unabsorbed flux in the 0.3-lOkeV band 
of ~ 3.5 X 10^* ergs ' cm ^ from ROSAT-HRI 1990 July 25 observation (Primini, Forman & Jones [T993t . ~ 1.0 x 10"'^ 
erg s"' cmr^ from Chandra-AClS 1999-2000 observations (Kong et fl/.'^UU^t, ~ 8.5 x lO"'"* erg s"' cm"^ from Chandra-HRC 
1999-2001 observations (Williams et al. 2003 1, ~ 3.4 x lO"'-* ergs"' cm"^ from Chandra-HRC 2001 October 31 observation 
(Kaaret I2002> . Note that both Primini, Forman & Jones J1993> and Kaaret (I2002> use a thermal bremsstrahlung distribution 
to fit the spectra, with temperature kT = 5keV and 2keV and absorption column A^h = 7 x 10^" cm"^ and 6.66 x lO-^cm"^ 
respectively, while Kong et al. (,2002 1 and Williams et al. (2003 1 use a power law with photon index F = 1.7 and absorption 
column A^H = 10^' cm"^. These values are within a factor five from the observed flux of ~ 1 .7 x 10"'^ erg s"' cm"^ we obtain from 
XMM-A^ewfon observations in the 0.3-lOkeV band (see Tableland sectionlT^. This indicates that XMMU J004308.6+41 1247 
is a fairly stable X-ray source, and it is not an X-ray transient. 



Mangano, V., Israel, G.L., Stella, L.: A possible new dipping X-ray source in M31 11 




Energy (keV) 

Fig. 7. phase-averaged spectra of XMMU J004308.6+41 1247 extracted from PN, MOS 1 and MOS2 data of the XMM-Newton 
observation carried out on 2002 January 6. The resulting simultaneous fit with the absorbed power law model is shown (see table 
O. The dashed line gives the best fit model for the PN spectrum, the dot-dashed is for the MOS 1 and the solid line for the MOS2. 




Energy (keV) 

Fig. 8. out-of-dip, dip and phase-averaged spectra of XMMU J004308. 6-1-41 1247 extracted from PN data of the XMM-Newton 
observation carried out on 2002 January 6. The resulting simultaneous best fit with the absorbed power law model with photon 
index F = 0.78 is shown (see Tables|3]and|3}. The dashed line gives the best fit model for the out-of-dip spectrum, the dot-dashed 
is for the dip spectrum and the solid line for the phase-averaged spectrum. 



5. Discussion 

The unprecedent throughput of XMM-Newton led to the detection of a 107 min modulation from the X-ray source 
XMMU J004308.6-H411247in the field of M31. This source is very likely the ROSAT-URl source [PFJ93]74 (Primini, Forman 
& Jones 1993 1. The modulation is confirmed also by Chandra archival observations of the source. The source flux was at a nearly 
steady level of ~ 1.7x10"'^ erg s"' cm"^ in the 0.3-lOkeV band throughout all XMM-Afewfon observations. This value is within a 
factor five from flux measurements obtained by other missions spanning the previous decade. The flux above corresponds to a lu- 
minosity of ~ 1 X 10"(ii/780 A;;?c)^ ergs where (/ is the distanceof XMMU J004308.6-I-41 1247 from us. The source spectrum 
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is well fit by an absorbed power law model with a photon index F ~ 0.8, but is also compatible with an absorbed Comptonization 
model {kTg ~ 3 keV and r ~ 13). No other simple single-component spectral model succeeds in fitting the XMM-Newton data. 
Owing to poor statistics the A^h value could not be measured from the XMM-Newton spectra and only upper limits were derived. 

The 107min modulation is due to a dip in the light curve that covers about one fourth of the cycle with an energy- 
independent intensity decrease compatible with 100% at dip minimum. The apparent smoothness of XMMU J004308. 6-1-41 1247 
light curves might be due to poor statistics. An eclipse-like behaviour is also compatible with the folded light curve of 
XMMU J004308. 6-1-41 1247. In this case, however, the "eclipse" duration would be of 0. 12 phase units during the XMM-A^ewfon 
observation of 2001 June 29, inconsistent with the ~ 0.2 duration obtained for the other three observations. This indicates that 
the observed modulation might not be compatible with a true eclipse by a companion star. 

The possibility that XMMU J004308. 6-1-41 1247 is foreground source in the Galaxy cannot be ruled out at present. Indeed 
the somewhat lower value of A^h, compared to the galactic A^h, that we inferred from the analysis of the X-ray spectrum would 
appear to favor this possibility. 

5.1. Galactic scenario 

For an assumed distance of 1 kpc, roughly coiTesponding to the thickness of the Galaxy in the direction of M31, the 0.3-lOkeV 
luminosity of XMMU J004308.6-)-41 1247 would be ~ 2 x 10^' erg s"', close to the low-end of (but still compatible with) the 
quiescent X-ray luminosity of magnetic cataclysmic variables (CVs; Ramsay et al. 1994 Ramsay & Cropper 1994 1. The period 
of the X-ray modulation of XMMU J004308.6-)-41 1247 is suggestive of the orbital period of a polar or intermediate polar (IP). 
Orbital periods of polars and IPs are mostly concentrated in the in the 1-10 hr range (Downes et al. 2001 1. CuiTent limits on the 
magnitude of any optical counterpart to XMMU J004308.6-H41 1247 (V > 21.7, R > 20.2, see Kong et al.Q^Ql. Haiman et al. 
I1994> are also compatible with a magnetic CV interpretation (see Downes et a/. l2001ii . 

Dip-like X-ray modulations and/or on/off shaped eclipses reflecting the orbital motion of the binary are seen in both classes 
of polars and IPs. Examples of this are Fo Aqr, AO Psc (Hellier, Garlick & Mason . 1993 ) and XY Arietis (showing eclipses, 
HeUier n"997> among IPs; 2A031 1-227 (EF Eri, Beuermann, Thomas & Pietsch n'991> . DP Leo (showing ecUpses, Schwope et al. 
120021 and, to lesser extent, AN UMa (Ramsay et aZ. 1 19941 1 among polars. 

However, the longest duration of an eclipse obtained through occultation of a compact accreting X-ray source by a donor star 
filling its Roche-Lobe of radius in a binary system with orbital period P, binary separation a and mass ratio q - Mjo„o,.IMx < 1 
can be calculated as j2ri/^ - ^ (in phase units). Using the relation of Eggleton (1983 1 to express the dependence of the 
Roche-Lobe radius of the donor on q we derive an upper limit to an eclipse duration of 0.4/7r ~ 0.13 phase units. A 0.2 
phase units eclipse of the X-ray source cannot be due to occultation by the companion star. Thus, in any interacting binary 
interpretation, the observed X-ray modulation of XMMU J004308 .6-1-41 1247 cannot be explained as due to a true eclipse. 

If the dip ingress/egress in the XMMU J004308.6H-41 1247 light curves sharp on the basis of the square-wave fit, a more 
plausible explanation for the observed modulation would be that of eclipses by the accretion column in a Polar. In this case no 
limit to the eclipse duration would be expected and occasional variations in the light curve profile may occur (see the case of EP 
Dra, Bridge etal. lM^ . 

The alternative scenario for XMMU J004308. 6-1-41 1247 X-ray modulation at the orbital period would be occultation of an IP 
by the plasma bulge due to impact of the accretion stream on the accretion disk. Such a mechanism, at work both in IPs and low 
mass X-ray binaries (LMXRBs), might produce a relatively long dip which, when averaged, has often a fairly smooth profile, 
consistent with the Gaussian fit to the folded light curves that we described in section BTTI The apparent energy-independence 
of the dips (to within the uncertainties) might be accounted for if the plasma cloud producing the dips had a very high degree 
of ionization, and the flux reduction were due to Thomson scattering rather than photo-electric absorption. In the IP scenario 
the ionization front due to the central X-ray source cannot reach the edge of the accretion disk: to produce the > 90% energy 
independent decrease of the source luminosity, an electron column Ne of more than 3.5 x 10^"* cm"^ would be required. In order 
to have this amount of ionized matter (with solar composition) in a cloud of thickness AR at a distance R from the central source 
emitting an X-ray luminosity Lx, the ionization parameter ^ = should be larger than 1000 (Tarter, Tucker & Salpeter 

119691 Hatchett, Buff & McCrav fl976> . We define e - i.e. the depth of the material that produces the dips expressed as a 
fraction of its distance from the source. We conservatively assume e - 1. Thus, in the case of XMMU J004308.6-(-41 1247 the 
distance from the central X-ray source of the material causing the dips should he R - < 10"^^. For an X-ray luminosity 
of Lx ~ 2 X 10^' erg s"' this maximum distance amounts to less than 1 km, i.e. orders of magnitude lower than the radius of any 
bulge in the accretion stream/disk interaction. 

Most IPs display soft X-ray modulation at the white dwarf (WD) spin period (33 s to 2 h, but always shorter than the orbital 
period, Patterson 1994 1 with energy-dependent amplitudes in the 10-40% range. No evidence was found for such a distinct X-ray 
modulation (possibly, in the soft X-ray range) resulting from the white dwarf rotation as observed in a number of IPs, in the case 
of XMMU J004308. 6-1-41 1247. We note, however, that X-ray pulsations with amplitude <30% would have remained undetected 
in the XMM-Newton fight curves of XMMU J004308.6H-41 1247. 
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A final possibility, producing a relatively smooth modulation in the X-ray light curve, is a self-eclipse of the hot emitting 
region on the polar cap of a highly magnetic WD (polar) in synchronous rotation with a low mass companion from which it 
accretes matter via Roche-Lobe overflow without accretion disk formation. In this case the profile and duration of the eclipse 
would be due to the shape and extension of the hot polar cap periodically hidden at the WD spin period, equal to the orbital 
period. The energy independence of the eclipse depth would be also explained. 

If XMMU J004308.6+41 1247 were a Galactic CV as suggested by its low X-ray luminosity for a 1 kpc distance scale, the 
folded light curve modulation would be clearly more suggestive of an accretion stream eclipse or a self eclipse of the WD in 
binary system of the polar class. 

However, evidence against the magnetic C V interpretation derives from the analogy with the spectral characteristics observed 
from Polars and IPs. In both cases the observed spectra are compatible with a thermal bremsstrahlung spectrum with temperature 
in excess of tens of keV, plus (in many cases) a very soft spectral component consistent with a blackbody of temperature in the 
10-50eV range. The bolometric luminosity in this component is from one to several hundreds times higher than the thermal 
bremsstrahlung component in Polars (Ramsay & Cropper 1994"i, while it ranges from zero (i.e. the component is not detected) to 
several tens of times the thermal bremsstrahlung component in IPs. 

The spectral characteristics of XMMU J004308.6+41 1247 are at variance with those of magnetic CVs. The X-ray spectrum 
of XMMU J004308.6+41 1247 has a characteristic temperature 2-3 keV (that of the Comptonising electrons) and, in any case, 
cannot be fit with a thermal bremsstrahlung spectrum. The addition of a soft blackbody component with temperature of 50 eV or 
less produces a bolometric blackbody flux comparable to the bremsstrahlung, but the fit should be discarded because of reduced 

of 1.6 (in the best case). 

Finally we note that, in the Galactic interpretation, the X-ray luminosity and spectrum of XMMU J004308.6+41 1247 would 
be similar to those observed from quiescent soft X-ray transient low mass binaries hosting a neutron star (see e.g. Campana 
& Stella EUnO)- X-ray eclipses have so far been observed only in a very low luminosity state of MXB 1659-29 (Wijnands et 
al. I2003t : no evidence for a dip-like modulation of the X-ray flux at the orbital period has yet been found in any quiescent 
neutron star soft X-ray transient. Therefore, while in principle viable, the possibility that XMMU J004308.6-(-41 1247 is a dip- 
ping quiescent soft X-ray transient is not especially appealing. On the other hand the interpretation of he X-ray modulation of 
XMMU J004308.6+41 1247 in terms of eclipses by the companion star is untenable because of the eclipse duration exceeds the 
maximum allowed duration (see above). 

5.2. Extragalactic scenario 

As anM31 source, the nearly steady luminosity of XMMU J004308.6-I-41 1247 (~ 10"ergs"' (0.3-lOkeV) would be similar to 
that of moderately bright LMXRBs (White, Stella & Parmar 1988). The shape of the 107 min periodic modulation, comprising 
a fairly flat on-phase lasting for about 3/4 of the cycle and a relatively long and smooth dip (1/4 of the cycle) reaching (nearly) 
zero flux over some 1/10 of the cycle, is clearly reminiscent of X-ray dipping LMXRBs seen from high inclinations. 

The possibility that the observed modulation is due to true eclipses of a compact central X-ray source by the companion star 
(supported by the square-wave fit of the folded XMMU J004308.6+41 1247 light curves shown at the end of section lTjI i. can be 
simply ruled out because the observed eclipse duration is incompatible with the geometry of the binary system (see section lOt . 

LMXRB dips are believed to arise from absorption and/or Thomson scattering of X-rays coming from a central source by a 
bulge in the outer regions of the accretion disk, where the stream from the Roche-Lobe filling companion impacts the outer disk 
rim. The dipping phenomenon reflects the orbital period of the system, as demonstrated by a few systems which display both 
X-ray dips and eclipses from the companion stars (e.g. MXB 1659-298, Cominsky & Wood 1984 Cominsky & Wood [T988l 
Wijnands et a/. [20031 EXO 0748-676, Parmar et a/. 119861 GRS 1747-312, into Zand et al. 2000 in't Zand et al. 2003 1. Despite 
some jitter in phase and the common presence of pronounced dip variability, the position of the bulge must remain (nearly) fixed 
in phase (Frank, King & Lasota 1987 1. It is estimated that the dip phenomenon occurs in systems with an inclination greater 
than ~ 70°. Some 20 dipping LMXRBs are currently known in the Galaxy. Only one of these (X 1916-053) has an orbital period 
shorter than ~2hr (White & Swank [T982t . while two (4U 1755-338 and 4U 1630-47, in Mason et al. 1985 and Tomsik, Lapshow 
& Kaaret|l998| respectively) likely host an accreting Black Hole Candidate (BHC). Most dipping LMXRBs are type I X-ray 
bursters and thus contain an accreting neutron star. 

In the dipping LMXRB interpretation XMMU J004308. 6-1-41 1247 would be a compact binary, the second dipper with a 
period below the period gap. 

The spectra of high luminosity (> 10^^ erg s ') LMXRBs hosting a central neutron star are usually well fit e.g. by a 
Comptonization model plus a blackbody component (White, Stella & Parmar 1988' for a review of other models see e.g. Di 
Salvo & Stella 2002). When required, the blackbody component has a temperature around 1.5 keV and it accounts for ~ 10-30% 
of the total 1-lOkeV flux. For both low/high luminosity sources, Comptonization parameters are y ~ 2-4, optical depths 
T ~ 10 - 15 and electron temperatures ~ 2 - 4keV. Spectra of (soft state) Black Hole Candidates in LMXRBs are well fit by 
Comptonization model too, but with a somewhat higher scattering depth (t > 15) and lower electron temperature (Tg ~1 keV; 
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White, Stella & Parmar'1988'). LMXRBs with somewhat lower luminosity (< 10^^ erg s ') are often dominated by a power law 
with photon index Til (Christian & Swank [T997l Schulz [T999l Church & Balucihska-Church 2001 1. 

In the case of XMMU J004308.6+41 1247, the XMM-Newton spectra ai-e reasonably well fit by a power law or a 
Comptonization model (see TableO. The power-law fit gives a photon index (F < 1), i.e. harder than that found in LMXRBs. 
The addition of a blackbody component to the power law, though not required from the statistical point of view, implies too high 
a blackbody temperature (~ 2.3 keV) and fraction of the total luminosity (~ 70%) as compared to typical values for LMXRBs. 
On the contrary the thermal Comptonization model discussed in section vields an electron temperature of ~ 3keV and an 
optical depth ~ 13 consistent with high luminosity LMXRBs hosting a neutron star. 

In this interpretation other relevant features of XMMU J004308.6+41 1247 are 

(i) the ~25% duty cycle of the dips; this is in agreement with those observed in Galactic dipping sources: the duty cycle of 
the dips varies from source to source and is typically 10%-40%. Sources showing dips with a ~20% duty cycle are, for example, 
4U 1755-338 (Mason et a/. fT985 ), XB 1254-690 (Courvoisier et a/. [T986t . MXB 1659-29 (Cominsky & Wood 1984i 

(ii) the shape of the dips, consistent with a relatively smooth average profile; this is common in Galactic dipping sources as 
well, despite the fact that individual dips often display pronounced irregular variability. Owing to poor statistics, individual dip 
variability would not be revealed in XMMU J004308.6+41 1247. 

(Hi) the very high degree of modulation, nearly 100%; in turn this is not unusual in LMXRBs: e.g. 4U 1916-053 and XB 1254- 
690 often display more > 90% modulated dips; MXB 1658-298 presented ~9Q% modulated dips in the Beppo-SAX observation 
of August 2000 (Oosterbroek et al. [2000l l: the dips of EXO 0748-676 occasionally reach a depth of 80% (Parmar et al. 198^. 

(iv) the absence of a conspicuous increase in A^h during dips. The most constraining upper limit to the average absorption 
column of the obscuring bulge producing the dips we got from XMM-Newton observations of XMMU J004308. 6+41 1247 is 
rather low: < 2 x 10^' cm"^ (at 95% confidence level; see Tablel^. However, we reiterate that the dip spectra were extracted over 
a fairly large phase interval (1/4 of a cycle), which samples predominantly the relatively shallow stages of the dip (this was done 
in order to increase the S/N). For this reason the upper limit we derived might underestimate the actual limit on the increase of 
A^H around dip minimum. 

(v) the energy independence of the modulation (to within the uncertainties). We have already discussed this point in section 
I5.1l (5th paragraph). As in the IP case, this may be accounted for if the disk bulge producing the dips were ionized by the central 
X-ray source and the flux reduction mechanism at work was Thomson scattering and not (as often happens) photo-electric 
absorption. 

Following the same line of reasoning of section 15m we conclude that for an X-ray luminosity Lx ~ I x 10^^ erg s ' the 
distance from the central X-ray source of the material causing the dips should be /? i 3 x 10^ cm, since only material closer than 
3x10^ cm from the X-ray source will be significantly ionized. By contrast, for a 107 min binary period the binary separation of 
the system is expected to be a ~ 5. 2x10^^ M^J^ (I +qy^^ cm, with Mx the mass of the compact object in solar units, and q the mass 
ratio of the donor star and the accreting star Thus, the accretion disk is expected to have a radius of ~ 0.3a i 2 x 10'"m],^^ cm 
for stellar masses with ratio q ^ 1 (Whitehurst & King '1991* Eggleton 1983). This requires, for ionization to be important, that 
the disk is at least one order of magnitude smaller than expected or the material responsible for the dips is located further inside 
the edge of the disk, both of which seem unlikely. We note that also in the case of 4U 1755-338 the ionization hypothesis leads 
to the same difficulties (Mason et a/. ll985l l. 

Another possibility to interpret the energy independent dips of XMMU J004308. 6+41 1247 would be the eclipse by the 
companion star of an extended source of X-rays of comparable size, such as that originating from electron scattering in an 
Accretion Disk Corona (ADC). Besides the energy-independence, this would explain the smooth ingress and egress of the dips 
as well as the (nearly) 100% reduction of the observed X-ray flux. In this hypothesis, however the central X-ray source should 
be completely hidden by the accretion disk rim in order to explain the absence of sharp eclipses. In analogy with known ADC 
sources we would thus expect that the true luminosity of XMMU J004308.6+41 1247 is one or two orders of magnitude higher 
than that observed, i.e. 10^"*-10^^ erg s"' . This in turn would suggest a BHC accreting at very high rates; however the observed 
source spectrum is considerably harder than that expected from such a BHC. 

6. Conclusion 

We discovered 107 min periodic dips in the X-ray light curve of XMMU J004308. 6+41 1247. The dip profile and amplitude, 
and average spectra of XMMU J004308. 6+41 1247 are consistent with a accreting neutron star in a short period LMXRB seen 
at a high inclination and located in M31. The apparent energy-independence of the dips is more difficult to interpret in this 
scenario. If this interpretation was correct, XMMU J004308. 6+41 1247 would be the second dipping source discovered in M31 
(after XMMU J0043 14. 1 +410724 in the globular cluster Bo 158, Trudolyubov et al. 2002 ). 

The possibility that XMMU J004308.6+41 1247 is a foreground magnetic CV cannot be ruled out at present, but appears to 
be less likely. 
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